Filtering random matrices: The effect of imperfect channel control in 

multiple-scattering 
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We present an analytic random matrix theory for the effect of imperfect channel control on the 
measured statistical properties of the scattering matrices of a disordered multiple-scattering medium. 
When the fraction of the controlled input channels, mi, and output channels, m2, is decreased from 
unity, the density of the transmission eigenvalues is shown to evolve from the bimodal distribution 
describing coherent diffusion, to the distribution characteristic of uncorrelated Gaussian random 
matrices, with a rapid loss of access to the open eigenchannels. The loss of correlation is also 
reffected in an increase in the information capacity per channel of the medium. Our results have 
strong implications for optical and microwave experiments on diffusive scattering media. 

PACS numbers: 42.25.Bs, 05.60.Cd, 02.10.Yn 



Describing wave propagation in strongly scattering 
media is a fundamental challenge in disordered sys- 
tems theory, relevant to problems in electromagnetism, 
acoustics and electron transport. The scattering ma- 
trix, which encodes fully multiple scattering within the 
medium, is a powerful tool, that relates arbitrary inputs 
to their outputs, and in principle allows the reconstruc- 
tion/prediction of either. For many years it has been 
understood that elastic multiple scattering can lead to 
Anderson localization of waves, and that even in the dif- 
fusive regime it creates important correlations in the scat- 
tering matrix of a random medium [l[ . In electron trans- 
port, the input electron state is uncontrolled and only a 
few statistical properties of the scattering matrix can be 
measured, through, e.g., conductance and shot noise ex- 
periments, but in classical wave systems it is possible to 
prepare specific input states and measure approximately 
the full S-matrix. There has been a great deal of interest 
in doing this in optical systems recently, since this knowl- 
edge allows the synthesis of input states, using spatial 
light modulators [2|, which provide dramatic control of 
the transmitted and reflected waves. Most strikingly one 
can focus light within or at the output of a strong scatter- 
ing medium [3|, |5| , with potential applications for imag- 
ing biological tissue, enhancing the sensitivity of spec- 
troscopy and, potentially, allowing the transmission of 
information through media which are opaque to typical 
input states [a lZJ . 

For most cases under experimental study, the S-matrix 
can be naturally divided into blocks containing trans- 
mission and reflection matrices, t and r, with a certain 
number of input and output channels, N . The possibil- 
ity of strongly enhanced transmission of waves through 
multiple scattering media was first discovered theoreti- 
cally more than twenty years ago [8 10], when it was 
shown that the eigenvalues, T„, of the Hermitian ma- 
trix tH have a bimodal distribution consisting of a large 
number of strongly reflected "closed" eigenchannels, and 
G "open" eigenchannels with T„ ~ 1 (where G — NijL 
is the dimensionless conductance, L is the sample length 



and € ^ L is the elastic mean free path) . If one were able 
to prepare the coherent superposition of channel states 
corresponding to an open eigenchannel, the input state 
would be transmitted with near unit efficiency through 
an effectively opaque medium (T = ^/L ^ 1). A distinct 
but related effect, discovered more recently, is coherent 
enhancement of absorption (CEA) [11|, the possibility of 
preparing an input state which is very strongly absorbed 
in a "white" medium with a £ <$^ L and, £a, the inelastic 
absorption length, greater than L. 

Despite these exciting predictions of theory, experi- 
mental efforts have as yet been unable to measure these 
effects. In particular, measurements of the transmission 
eigenvalue density have not revealed the predicted bi- 
modal structure. The eigenvalues of the Hermitian ma- 
trix tH are the square of the singular values, \/T^, of the 
complex matrix t, the distribution of which was first mea- 
sured in acoustics [12| and then in optics [J] , with results 
very close to the "quarter-circle law" |13|, characteristic 
of uncorrelated Gaussian random matrices. The discrep- 
ancy with the theoretical prediction was attributed to 
the incomplete angular coverage of the input and out- 
put channels in the experiments [2,]. Efforts have been 
made to measure a larger fraction of the transmission 
matrix in order to reveal the existence of open eigenchan- 
nels |14| . but strikingly, even in experiments performed 
with microwaves in a multimode waveguide, where al- 
most all channels can be addressed, the measured distri- 
bution does not reveal the second peak associated with 
completely open eigenchannels 15 1. 

These observations motivate the study of a new ran- 
dom matrix ensemble, directly relevant to the experimen- 
tal measurements, which we refer to as a filtered random 
matrix (FRM) ensemble. The definition of this ensemble 
and many of our results are quite general, and would ap- 
ply to S-matrices in arbitrary scattering geometries, and 
more generally, to any situation in which only a portion 
of a physical random matrix. A, is measurable; however 
here we apply the theory to the important case of ran- 
dom transmission and reflection matrices. In almost all 
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FIG. 1: Transmission eigenvalue density of a disordered slab 
placed in a waveguide with A^ = 485 channels (length L = 
150/A:, width W — 900/fc), for different fractions of controlled 
channels m — M/N . Numerical results (dots) are obtained 
from solving the wave equation for 120 realizations of the slab, 
with dielectric function e(r) = rig + 5e(r); no = 1.5 and 5e(r) 
is uniformly distributed between [—1.1, +1.1] with e^c = 1-1 
in the slab and 5e = Q m the empty waveguide. The solid 
lines are the theoretical prediction based on Eqs. ((2]), ((3| and 

(0, where T = (E^Li^")/^ = 0-067 is found from the 
simulation with perfect channel control (m = 1). Inset shows 
the maximal transmission enhancement possible for a given 
m, where (T) = mT, and (Tmax) is calculated by method 
given in [lq |. 



such experiments there are some hmitations on generat- 
ing input channel states (e.g. due to the limited numer- 
ical aperture associated with the incoming light) and in 
detecting the outgoing states; we will refer to both situ- 
ations as imperfect channel control (ICC). Hence the ex- 
periments will typically not have access to the full, NxN, 
matrices, t, r, but rather to some finite sub-matrix of t, r. 
We derive below the statistical properties of such mea- 
sured sub-matrices, to assess the effects of ICC on the 
correlations inherent in diffusive transmission and reflec- 
tion. 

The effect of ICC is to map the NxN matrix A to 
A = P2AP1 , where Pi , P2 are N x Afi and M2 x N matri- 
ces which eliminate N — Mi columns and N — M2 rows, 
respectively, of the original random matrix A. A is the 
measured random matrix with ICC and Mi,M2 < N 
are the number of input (output) channels controlled. 
We will compute the eigenvalue density p^',^{x) of the 
Hermitian filtered random matrix, A^ A^ in the limit 
N,Mi,M2 — )■ 00, with arbitrary but fixed fractions of 
controlled channels, mi = Mi/N and 1112 = AI2/N, as- 
suming that the eigenvalue density (and resolvent, see 
below) of the full matrix. A'' A, is known. 

In order to obtain the eigenvalue density p^^^{x), it is 
convenient to introduce the resolvent 



^-'-^^^ = M^'^^) 



where averaging (...) denotes the ensemble or disorder 
average. The eigenvalue density is given by 

PAiAi^) = lim Im5Ati(a^ + i*?)- (2) 

Although field-theoretic or diagrammatic approaches to 
this problem are possible, here we take advantage of the 
multiplicative structure of A^ A = P\Al^P\P2AP\ to em- 
ploy the powerful method of free probability theory [l6| . 
This theory identifies a sufficient condition — "asymp- 
totic freeness" — under which the spectral properties of 
a product of matrices can be found algebraically from 
the spectral properties of the factors. Loosely speaking, 
asymptotic freeness can be thought of as the general- 
ization of statistical independence to the case where the 
random "variables" do not commute |17|. In our case, A 
is assumed to be random, and the matrices P\ and P2 can 
be generated by randomly suppressing rows or columns of 
a,n N X N identity matrix. If the "lost" channels are not 
chosen randomly, the eigenvalue density is still given by 
the present theory with renormalized parameters, which 
can be determined from microscopic treatments of the 
radiative transfer equation [18]. Here we only discuss 
ICC in channel (momentum) space, but we find that the 
FRM distributions we obtain can describe the behavior 
of a focused spot of light (real space filtering), incident 
on a waveguide or even on a slab with no walls at all [13] . 
Specializing general results of free probability theory to 
this specific ensemble, we show in [^18] that the unknown 
resolvent g^^^{z) may be obtained from the known re- 
solvent gA'TA by nieans of the implicit equation: 

N{z)gA,A{Niz)yD{z))^D{z), (3) 

where N{z) and D{z) are two auxiliary functions defined 



N{z) = zmig^;^{z) + 1 - toi, 

D{z) = mig^i^{z) [zmig^^^{z) + ma 



mi[ 



(4) 
(5) 



We obtain the results given below by solving this self- 
consistent equation in the complex z-plane numerically 
and taking the limit of Eq. ([2]) . The theory also gives us 
explicit formulas for moments of the eigenvalue density 
of A^ A in terms of moments of A^ A, which we use at 
some points in deriving the results given below |18| . 

We now use Eq. ^, setting ^ = t, to study the effect 
of ICC on the transmission through a disordered non- 
absorbing slab in the diffusive regime, in which NT — 
G > 1, where T = (^^=1 Tn) /N = ijL ^. The resol- 
vent associated to the bimodal transmission eigenvalue 
"113, [3 is 



density p^t* 
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zArctanh 
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(1) 



The solution for the density p^f^{T), obtained from 
Eqs. (|3]), ([S]) and ([2]), is shown in Fig.[T] where we chose 



mi = 1712 = TO. For m = 1, Pifi{z) has the expected 
bimodal shape, even for the slab geometry simulated 
in Fig. m confirming that this distribution is not re- 
stricted to the quasi-one-dimensional geometry L ^ W 
[20| . The number of open eigenchannels (the channels 
with T > 1/e) is equal to the dimensionless conductance 
G = NT; for the most open eigenchannel, (T™""^) -^ 1 
as N ^ oo. Introducing a small degree of ICC (m < 1) 
abruptly suppresses the most open eigenchannels: the 
mean of the largest eigenvalue l^T"^""^) becomes strictly 
smaller than 1 as A'^ — >■ c», and the distribution loses its 
second characteristic peak. This striking property (pre- 
served for TOi 7^ 777.2) indicates that phenomena based 
on extremely open eigenchannels are highly sensitive to 
ICC and it may explain why the bimodal shape has not 
been observed in real experiments, even with almost per- 
fect channel control [15|. Note, however, that even if the 
bimodal shape is lost, (T™"^) > T can still hold (sec 
inset to Fig. [ij, for reasonable values of 7r7, so strongly 
enhanced total transmission is not ruled out by ICC. 

Our analytical prediction is in excellent agreement 
with the result of numerical simulations of the wave equa- 
tion [V^ + fc^e(r)]7/)(r) = 0, based on numerical discretiza- 
tion in a two-dimensional disordered slab embedded in a 
multimode waveguide with iV = 485 channels and per- 
fectly reflecting walls. The N x N transmission matrix 
is computed using the recursive Green's function method 
2l| , and members of the filtered ensemble are then gen- 
erated by random projection. 

When m is further reduced, the correlations con- 
tained in the transmission matrix are progressively lost 
and pm{T) evolves such that the distribution ol X ~ 
yjT I (T) converges to the quarter circle law, pjc(x) = 
v4^-a?/7r, independent of i. Thus universal, uncorre- 
lated behavior is reached when 777 < T (or M < G), 
in agreement with measurements reported in [J]. This 
loss of correlations in the limit of small degree of channel 
control remains when 7711 ^ 7772, but in a more subtle 
form. For example, for 7771 ^ m < T and 7772 = 1, 
we find that px approaches the Marchenko-Pastur (MP) 
law [l3|, describing rectangular random matrices with 
uncorrelated Gaussian matrix elements, but for a ma- 
trix ensemble with a disorder-dependent, effective value 
of 777 — >■ 777. Specifically 



100 



Px{x) 



1 



■■\/(.T+ — x^) (a;^ — a; ), 



(7) 



where x 



1± 



and m = 777(2/3r- 1) [191; this 

corresponds to a MP distribution for N x M matrices, 
withiV = 3TA/'/(2-3f'). 

The statistics of lossless reflection with ICC can be 
obtained similarly to those of transmission, with qualita- 
tively similar results (i.e. suppression of extremal val- 
ues and convergence to an effective MP distribution). 
However, in the case of an absorbing disordered medium 
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FIG. 2: Reflection eigenvalue density of a disordered ab- 
sorbing slab with the same geometry and dielectric function 
as in Fig. [l] except for the addition to e(r) of a constant 
imaginary part, 0.0031, representing absorption. A fraction 
m = M/N of the input channels are excited while all output 
reflection channels are collected. Numerical results are based 
on 100 realizations of the disordered slab. Solid lines are the 
theoretical prediction based on Eqs. ([2]), ((3]) and ((9|; where 
a = 0.023 is determined by the numerical value of _R = 0.74 
with m = 1. The ballistic and diffusive absorption lengths 
are la = £/a = 2.95L and ^/IL = 0.44L. 



within a waveguide, one has a distinct statistical ensem- 
ble 22, l23[ from that in the lossless case. The extremal 
eigenvalue statistics of this ensemble were recently stud- 
ied by Chong and Stone [ll[ and lead to the phenomenon 
of coherently enhanced absorption (GEA). Here the (non- 
unitary) S-matrix and the reflection matrix coincide, and 
1 — Rn represents the absorbed fraction of an incident 



eigenchannel, where {Rn} are the eigenvalues of rV. Let 
^a ^ ^ be the ballistic absorption length and consider 
the regime in which i <^ -JWa < L,£a, so that elastic 
scattering is strong, absorption is weak, but transmis- 
sion is negligible. In Ref. [llj it was found that when 
N'^{i/ia) ^ 1, the smallest Rn (reflectivity of the most 
highly absorbed eigenchannel) was orders of magnitude 
smaller than the mean reflectivity, R: 



(iJ™") 



1 



R 



2N^c 



-«^- 

<■(!, 



(8) 



As A^ ^- 00, (i?™™) — >• while R remains '^ 1, which 
would not be true, e.g. for the MP law, and is the essence 
of CEA; in addition the density p^tr(-R) diverges at i? = 
(see Fig.©. The l/N"^ scaling of (i?"*™) holds even when 
the absorption is non-uniform, e.g. for a buried absorber 
behind an "opaque", lossless layer. The effect of ICC 
in this case is again found by solving Eq. ([3]), now with 
A — r. We will specialize to the case where a fraction 
rr7 — M/N of the input channels can be excited, while 
the field in all output channels is collected, 7771 = tr, 
m2 = 1. We find the eigenvalue density Pf\f^ from the 



known density p^t^ 23, l23| and the associated resolvent, 



z - 1 + 2a - 2a J\ + 1/a - Xlaz 

9rir{z) = {1-z)^ ■ ^^ 

Our results (Fig. [5]) again show excellent agreement with 
the simulation of the relevant wave equation for any m. 
Of particular interest is the behavior and support of the 
density near i? = 0. If we take the limit A^ — >■ oo and 
then consider ni — 1 — (5, we find that the density has no 
support at i? = 0, but instead the support has a sharp 
cut-off at i? ^ 0, which we can identify with (i?"""). 



A general equation for (i?™™) is derived in [18|, and 
its solution, normalized by (R) = R, is plotted in the 
inset to Fig. H The strong sensitivity of (i?™™) to TV 
is completely lost for all m such that S > 2^/2/N, and 
instead we find to leading order, 



(5^ 
^ ^ 16a 



(10) 



for (5 < VS, and (i?""') 



6 > 



7,1/3 



(5 - 3(1 - ,5)i/3,52/3ai/3 for 
The experimentally observable decrease in 
reflectivity relative to R will in most cases not be de- 
termined by N, but instead will be controlled by 6 and 
typically will be much less than predicted by Eq. ([8]). 
Note, however, that it can still be substantial, e.g. a fac- 
tor of ^ 5 decrease in reflectivity, when only half of the 
channels are controlled for realistic experimental param- 
eters. 

The joint probability distribution corresponding to this 
system, for m = 1, is that of the Laguerre ensemble of 
random matrix theory llj. The dramatic change in the 
support of the distribution p{R) at _R = from infinite 
to zero when m < 1 is related to a phase transition in 
the Dyson gas model as m is decreased from unity [l8J. 
A similar transition happens for p{T) near unity for the 
non-absorbing system. We observe in both cases that 
ICC changes the short-range correlations of the eigen- 
values, inducing a transition from linear (/3 — 1) eigen- 
value repulsion to quadratic (/? = 2) eigenvalue repulsion, 
similar to that normally associated with time-reversal 
(TR) symmetry breaking. This is due to the fact that 
randomly suppressing rows or columns of the complex 
matrices t, r leads to an ensemble of S-matrices which 
violate the usual TR symmetry constraint, SS* — 1. 
This suggests that in many experiments with nominal 
TR symmetry the T and R spectra will nonetheless 
show quadratic eigenvalue repulsion, and there is some 
evidence to this effect 2J]. Finally, for tti < a, the 



density of normalized absorption of each eigenchannel, 
X = v/(l - i?)/(l - -R) is of the uncorrelated form ^ 
with m = m(l + 2a) / A^J a{l + a) - m/2 ^. 

The previous analyses suggests that the degree of cor- 
relations contained in the matrices i, f is controlled by the 
parameters m/T, m/ a, so that when the fractional con- 
trol is less than the "loss rates" T and a, correlations are 
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FIG. 3: Information capacity per channel of a disordered 
multimode waveguide for a signal-to-noise ratio at the output 
SNR = 100. Solid curves are obtained from Eqs. (fTTI) . Q, 
©and (111), for different values of m = M/N and f = £/L. 
Dots correspond to the capacity p2|) with perfect channel 
control (m = 1) and the horizontal dashed line denotes the 
capacity used in wireless communications (see text), C — 5.48 
for SNR = 100. 



lost. To make this statement precise in the sense of infor- 
mation theory, we study the information capacity, C, of 
a disordered multimode waveguide, focusing on the case 
of transmission without absorption. The information ca- 
pacity C is the maximal rate, expressed in bits per second 
per Hertz (bps/Hz), at which the sender can transfer in- 
formation with a vanishingly low probability of error [25| . 
Correlations reduce C, because they reduce the range 
of possible outputs as the inputs are varied. More pre- 
cisely, C is the maximum of the mean mutual information 
carried between the input and output [17J, |25|, |26[ . Let 
ip™ be a vector of size Mi containing some information 
encoded on the A/i controlled channels at the input of 
the guide. The received signal in the M2 output chan- 
nels is V'°"' = t'i/''" + V'"- The vector ■0" accounts for 
an external Gaussian white noise: E[^^ip0*] = AfSap, 
where £^[. • ■ ] denotes the averaging with respect to the 
probability p(V''") that the sender transmits i/)™. For 
technical reasons we specialize to the case of a long dis- 
ordered waveguide, as opposed to a wide disordered slab, 
but we expect our conclusions to by qualitatively valid 
more generally. In this case, assuming uniform average 
incident flux in each channel llOl l. C can be written as 

M 



/■ 



C = Ml / dxKtt7(T)(«)log2(l + SNRx). 



(11) 



In this expression, SNR is the signal to noise ratio at 
the_output: SNR = P^'/A/", with P°"* = (T) P"' = 
77^2TP™. Our microscopic theory for Piif{x) allows us to 
compute C for arbitrary choice of mi and m^ [18| | . 

In Fig. [21 we plot the capacity per channel C /Mi for a 
given value SNR = 100 and mi = m,2 = m. Two limits 




can be identified. First, when one has fuU channel control 
(m = 1), we find 



(12) 



showing that (up to logarithmic corrections) the number 
of open eigenchannels G can be interpreted as the bi- 
trate of the disordered sample with perfect channel con- 
trol. Even though one measures N output channels, their 
correlations reduce their information capacity to G ^ N . 
Second, in the regime of strong ICC, m <T, C increases, 
becoming independent of T, and is given by C/M = 
log2 (1 + 4 SNR) - 2 - l/ln2 + (V1 + 4SNR - l) /21n2 SNR. 
This is the standard form used to model free space com- 
munication where many channels are uncontrolled |17| . 
In the intermediate regime, f'<TO<l,f'<0.1, C/M 
depends only on the ratio M/G, confirming that this ra- 
tio is a measure of the degree of correlations. As long 
as M > G the disorder-induced correlations are revealed 
and the capacity per channel drops from its maximum, 
uncorrelated value. 

In summary, the extremal eigenvalue properties neces- 
sary for transmission through opaque media or enhanced 
absorption are suppressed substantially as the degree of 
channel control is reduced, however strong enhancements 
should still be possible for achievable values of mi,m2- In 
most cases, experiments will need to measure mi and TO2 
in order to estimate the maximum enhancements possible 
for a given system and set-up. 
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SINGULAR VALUE DENSITY OF A FILTERED 
RANDOM MATRIX 

Let us consider a. N x N random matrix A such that 
the eigenvalue density PA'tAi^) oi the Hermitian matrix 
A' A is known. Since the singular values of A are the 
square root of the eigenvalues A of A'' A, the singular 
value density of A is p^ {x) = 2xpa'ia{^'^)- Our goal is 
to show how the singular spectrum of A is modified when 
we suppress a certain fraction of rows and columns in the 
latter. 

The modified spectrum is the one of a M2 x Mi sub- 
matrix of A (Ml , M2 < N) that can be represented as 
A — P2AP1, where Pi and P2 are two rectangular ma- 
trices that select Mi columns and M2 rows of A, re- 
spectively. In order to find the eigenvalue density of 
A'' A = pIA''P^P2APi, we introduce the resolvent 



gAiAiz) = 



Tr- 



1 



A'fA 



The eigenvalue density oiA'^A, 



Pa^a( 



1 / "" 



\n=l 



is related to the resolvent through 



Pa^a{x) = -- lim lmgAiA(x + iv), 

TT r(-i.O+ 



or, inversely, 



9AiA{z) 



dx 



PA^Aix) 



(1) 



(2) 



(3) 



(4) 



The combination of the definition ([TJ with the relation 
Tr(z - BC)-^ = Tr(z - CJ5)-i + {M - N)/z (valid for 
arbitrary matrices B and C of size M x N and N x M, 
respectively) yields 



9AfAiz) = 



9Air2AVi{z) I -mi 



mi 



miz 



(5) 



where toi = Mi/N, and Vi = PiP^ and V2 = P2P2 
are two projectors of size N x N. The problem is then 
reduced to finding the resolvent of the matrix A^'P2A'Pi. 
As explained in the main text (MT) , our idea is to use 
free probability theory [l|, |2[ . The eigenvalue distribution 
of the product of two asymptotically free matrices Xi and 
X2 can be found using the so-called 5-transform [2|. If 
we define x{z) as a solution of 



1 



X{z)' \x{z) 



l^z, 



(6) 



then the iS-transform is 

Siz) 



l + z 



X{z). 



(7) 



The iS-transform of the product X1X2 satisfies [2|: 

SxiX2{z)^Sxt{z)SxoAz)- (8) 

Since the projectors Vi and V2 are not random, we will 
assume that they are asymptotically free with respect to 
the random matrix A'^ A, in the limit N,Mi,M2 — > 00 
with mi = Mi/N and m,2 — M2/N fixed. Similarly, 
we assume that A^'P2A and Vi are asymptotically free. 
With the property ([8]), we obtain 



'Sa->V2AVi{z) = S-p2{z)SAiA{z)S-pi{z), 



(9) 



where we used SAiv2A{z) ~ ^^At-pgl^) ^^^ SAA^iz) — 
Saia{z)- In addition, the projector Vi {V2) has Mi (M2) 
eigenvalues equal to one and N—Mi {N—M2) eigenvalues 
equal to zero. Therefore, g-p^ 3(2) = m,i,2/{z — 1) + (1 — 
mi^2)/z. According to the definitions ([6]) and ^, the 
iS-transform of the projectors are: 

SvU^) = ^-^^- (10) 

z + mi^2 

Finally, combining Eqs. ^ and (ITUl) . and after some al- 
gebra using Eqs. ([S]), © and (O, we obtain Eq. (3) of 
the MT, which is an implicit equation obeyed by the un- 
known resolvent g^^^{z). If the singular value density of 
A is known, the singular value density of A is obtained by 
applying the following recipe. (1) Compute the resolvent 
of A^A with Eq. ([4]). (2) Find (numerically) the solution 
g^^^(z) of Eq. (3) of the MT. (3) Deduce the singular 
value density p^ (x) = 2xp^\^{x^) from Eq. ([3]). 

The different properties of the solution of the equa- 
tion obeyed by g^,;^{z) will be detailed elsewhere [3^. 
Here, we simply mention some properties implicitly used 
throughout the Letter: 

• Since A is of size M2 x Mi , the Mi x Mi matrix A'^ A 
is of rank min(Mi, M2), so that A^ A has max(Mi — 
M2, 0) eigenvalues equal to zero. 

• The eigenvalue density of AA^ is simply related to 
the one of A'^ A by 



PAAt(2^) = —PA^Ai^) 
777.2 



1 - ^ ) K^y 

"72/ 



(11) 



• The first two cumulants of the distribution Pa:\ A 
are 

(A^ti) = ™2 (A^-u) , (12) 

7772VarA^tA + (1 - ^2) (A^t^)^ 



VarA JtA = mim2 



(13) 
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FIG. 1: Variation of the smallest reflection eigenvalue 
(i?"'") with mi = m = M/N and a = t/ta- Each data point 
is an average over 100 realizations, with the same parameters 
as in Fig. 2 of the MT; a is determined by the numerical value 
oiR= (En=i -R") 1^ = l + 2a-2\/a(l + a), and solid lines 
follow from the equations for the edge of the distribution p^tf 
(see text). 



• It can also be shown that the edges x* of the dis- 
tribution are given by 



a; = 4 



nil — 1 



e<?AtA(e) 



m2 - 1 

e5AtA(e 



where ^* is the solution of 



d5AtA(0 



d^ 



5AtA(r 



(14) 



(15) 



e 2e 

-(l-mi)(l-m2)+r'gAtA(r)^ 

(1 _ ™^)(i _ ^2) _ (1 _ ^^/2 - m2/2)e5AtA(e) ■ 

Equations (|14p and ([T5|) have been used to obtain, in 
the limit A^ — > 00, (T™°^) without absorption (see in- 
set to Fig. 1 of the MT) and (i?™™) with absorption 
(see inset to Fig. 2 of the MT). In order to illustrate 
the accuracy of those predictions in the context of coher- 
ent enhancement of absorption, (i?™™), evaluated from 
Eqs. dUl), (HSl) and Eq. (9) of the MT, is compared in 
Fig. 1 with the result of numerical simulations of the 
wave equation, for different values of absorption. 



FILTERING THE TRANSMISSION AND 

REFLECTION MATRICES WITH A 
NON-RANDOM SUBSET OF CHANNELS 

In a diffusive sample, waves lose the memory of their 
initial direction after a distance of the order of the trans- 
port mean free path (.tn meaning that, inside the sample, 
all momentum channels diffuse in the same way. This 
observation justifies the fact that the projectors Vi and 
V2 should be asymptotically free with respect to iH and 
f'l^f. However, since the different incoming channels are 
not converted into diffusons at the same depth inside the 
sample, boundary conditions make the channels slightly 



non-equivalent. To illustrate this point, let us consider a 
waveguide with N propagating channels. A typical ex- 
perimental constraint might be that all the channels in a 
given angular range (numerical aperture) are controlled 
on incidence and collected on output, so the "lost chan- 
nels" are not chosen randomly. Using the radiative trans- 
fer equation for the specific intensity of light, it can be 
shown that an incoming channel oriented with an angle 
Oa with respect to the longitudinal axis of the waveg- 
uide is transmitted in the outgoing direction dfj with an 
amplitude {Tab) = (|Up) of the form [li] 



{Tab 



XaXb 



^N 



T 



(16) 



where Xa = {Ta) jT - Eb=i (^afc) /T. In 2D, Xa = 
2[cos(6'a) + zo/^trl/TT, and in 3D, Xa. = 3[cos(6'a) -I- 
zo/^tr]/4. Here, zq is the extrapolation length, equal to 
7rftr/4 in 2D and 2^4^/3 in 3D {Hj. It then follows that 
the mean value of the distribution p^tf is given by 



Ml A/2 



(T) 



Y^'^^^'^b) ^^^2S{b]{'m2)s{a]{'mi)T, (17) 



Ml 



a=l &=1 



,M 



where Si^a}{''n) — X]a=i Xa/M satisfies s^a}{'^) = 1- If the 
subset {a} is chosen randomly among the N channels, 
we get S{a} ('Ti) = 1 for arbitrary m and (T) = m2T, 
in agreement with Eq. (|12|) . The latter property is not 
exactly satisfied for a non-random subset {a}. If we select 
in 2D the Af channels with angles \9\ < Arcsin(TO), we 
would get s{m) = rn/2+Arcsm{m)/Tr+m\'l^^rn^/Tr > 1 
and (T) > mzf". 

The effect of the surface layers where free propagation 
is converted into diffusion can be taken into account by 
a renormalization of the ICC parameters. Equation (J17p 
suggests to renormalize the parameters TO2 and T as: 



7712 ^m2S{fc}(TO2), T -^ S[a}imi)T. 



(18) 



To test the validity of our theory, we have reevaluated the 
transmission and refiection eigenvalue distributions stud- 
ied in Figs. 1 and 2 of the MT with projectors correspond- 
ing to a limited numerical aperture \9\ < Arcsin(777), and 
compared them with ICC predictions. When we use the 
bare values of 7n2 and T, a shown in Fig. 2, one still finds 
reasonable agreement, but with small deviations from nu- 
merical data (solid lines). These discrepancies are nicely 
corrected by the renormalization (TTS)) (dashed lines) , con- 
firming that their microscopic origin is the surface layer 
effect studied above. 



PHASE TRANSITION IN THE DYSON GAS 

As explained in the MT, we observe a dramatic change 
in the support of the distribution pm{T) at T = 1 from 
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FIG. 2: Effect of a limited numerical aperture on the mea- 
sured transmission and reflection eigenvalue densities. Nu- 
merical results (dots) in (a) and (b) are obtained for the same 
parameters as in Figs. 1 and 2 of the MT, respectively. But 
here the projectors Pi and P2 are non-random and mimic a 
numerical aperture \6\ < Arcsin(m). (a) mi — m2 = m. 
Solid lines are identical to Fig. 1 of the MT, and dashed lined 
correspond to ICC theory with renormalized parameters (|18|l ; 
s{m) = X^a=i i'^'^) /MT is determined from the numerical t- 
matrix with m = 1. (b) mi = m, m2 = 1. No significant 
deviation is observed with respect to the results presented in 
Fig. 2 of the MT. 



infinite to zero when m < 1. A similar transition happens 
for pf\f.{R) at i? = 0, with or without absorption. To get 
more insight in this transition, let us first consider the 
case of an absorbing quasi-lD disordered waveguide. For 
TOi — 1712 — 1 and L > \/Ua^ the joint probability density 
(JPD)of the reflection eigenvalues i?„ G [0, 1] is of the 
form Q, 



where 



ui(i?„) = (iV + l) 



2a- 



1 - i?r, 



, «2(-R„,-Rm)] 



In (1 - i?„ 



U2{Rn,Rm) — — ln|i?„ — Rn 



(19) 

(20) 
(21) 



and /3 = I (unbroken time-reversal symmetry). The 
JPD coincides with the Gibbs distribution of a two- 



dimensionsionnal Coulomb gas with coordinates i?„ con- 
fined to the segment [0, 1] and experiencing the one-body 
potential ui. The Coulomb repulsion tends to spread the 
eigenvalues on the real line but because of the "hard wall" 
at i? = 0, they accumulate in the vicinity of the latter, 
giving rise to the peak of open eigenchannels observed 
near i? = in Fig. 2. On the other hand, for mi = m < a 
and TO2 = 1, the density oi X = \/(l — i?)/(l — R) is of 
the Marchenko-Pastur type [see Eq. (7) of the MT] with 
m = m(l-h2a)/4^a(l + a)-m/2 Q. The JPD of the M 
remaining eigenvalues corresponding to this law is still a 
Gibbs distribution, with [sl M 



Ui{Rn) 



M 
2^1 



1 — Rn 
1-R 



+ (m- l)ln(l-i?„) 



(22) 



1*2 unchanged, and /3 — 2 (broken time-reversal symme- 
try). In the intermediate regime a < ?n < 1, it is rea- 
sonable to think that the potential ui evolves between 
the two limits (l20t and (|22|) . so that the modification in 
(i?™™) can be interpreted as a zero-temperature phase 
transition in the associate Coulomb gas, occurring at 
m — 1 and due to a modification of the potential induced 
by imperfect channel control. 

A similar analysis can be performed for (T"^°-^) in a 
quasi-lD disordered waveguide. For mi ^ m,2 = 1, the 
JPD of the transmission eigenvalues has still the form of 
the Gibbs distribution dTOl) with 



Ul{Tn) 



N + 1 



[in T„ + Tx 



^1 



(23) 



U2{Tn, T„0 = -- [In |T„ - r,„| + In \xl - xU] , (24) 



Xn — tanh^ y/1 — Tn, and (3 = 1. The eigenvalue re- 
pulsion (|24p coincides with the Coulomb interaction (PT|) 
only for the most open eigenchannels r„ — > 1, while for 
the closed eigenchannels T„ — >■ 0, the interaction is twice 
as small |10| . When the fraction mi = to of input chan- 
nels is reduced, some eigenvalues are lost and the remain- 
ing eigenvalues redistribute on the real line. In the limit 
TO < r, we find that the distribution oi X = ^Tj (T) 
is given by Eq. (7) of the MT with m == m{2/'Sf - 1), 
suggesting that Pi has the same spectral properties as a 
complex Wishart matrix. The JPD of the latter is of the 
form (IT^ with 



Ul{Tn) 



M 
2^ 



Tn 

T 



+ (to — l)lnr„ 



U2{Tn, Tm) = - In |T„ - T„ 



(25) 
(26) 



and /3 = 2. The comparison of Eqs. ^^ and (|24l) with 
Eqs. (I25]) and (|26)) suggests that ICC induces a modi- 
fication of both ui and U2, in such a way that a phase 
transition — responsible for the abrupt change in (J"™"^) 
— occurs at 771 = 1 and the usual Coulomb interaction 
is restored for to < T. As noted in the MT, ICC is also 



characterized by a rapid transition /3 = 1 ^- /3 = 2, that 
has a clear signature in the unfolded level spacing dis- 
tribution: we observed numerically that the latter, given 
by the Wigner surmise typical of GOE for rrii = 7712 — 1, 
obeys the Wigner surmise for GUE in the case of ICC 
i- 



INFORMATION CAPACITY OF A DISORDERED 
MULTIMODE WAVEGUIDE 

Imagine a communication medium made of a disor- 
dered multimode waveguide with N propagating chan- 
nels, such that we control only Mi < N channels at the 
input and measure M2 < N at the output. As explained 
in the MT, the standard linear Gaussian model of such 
a system is 



^out ^ ^^^• 



r 



(27) 



where V'™ is a vector of size Mi that contains the complex 
amplitude of the input field expressed in the mode basis 
of the empty waveguide, ■i/'°"* is the output field of size 
M2, t is the M2 X Ml transmission matrix and ■0" is an 
external Gaussian white noise: -E[f/'2V'^*] = J^^ap- In 
addition, we assume that the input field obeys the power 
constraint 



TrQ s^ MiP* 



(28) 



where Q is the covariance matrix of the input, Qa(j = 

The information capacity of this medium is obtained 
by maximizing the mutual information T between ?/;"' 
and V'™* [11] ■ The latter is defined as 



/i(^««t)_/i(^°«t|^™)^ 



(29) 



where /i(V'°"*) = - /p(V'°"*)log2[p(V''"'*)]df/'°"* is the dif- 
ferential entropy of the output field and /i(-0°"*| ?/;*") — 
-/p('(/'™,V'™*)log2[p(V'°"*IV''")]dV''"dV'°"* is the condi- 
tional differential entropy. For a noise ■0" indepen- 
dent of the input, the mutual information reduces to 
I = /i(V'°"*) - h{tp"-). It can be shown that the in- 



put field that maximizes I must be Gaussian llj, so 
that the output itself is Gaussian, i.e. p{'ip°'^*) — 
det(7ri^~^)e~''' ^ "^ , with K the covariance ma- 
trix of 1/)°"*. Its entropy is /i(i/)°"*) — log2det(7reiir). 
Since -0" is also a Gaussian random field, we get X — 
log2det(if/A/'). The covariance matrix K can be evalu- 
ated from (P7|) . yielding 



X = log2det Im; 






(30) 



When the receiver at the output of the waveguide has 
only access to the statistics of t, the capacity averaged 



over the disorder configurations, C, is found by maximiz- 
ing the average mutual information over the distribution 
p{ip^^), i.e. over the covariance matrix Q y. 



C = 



(31) 



under the constraint (^5)1 . Using the identity det(/M2 + 
AB) = det(/Mi + BA), valid for arbitrary matrices A 
and B of size M2 x Mi and Mi x M2, we can express C 

as [il[il 



C — max ( log, dot 
Q 



'A/i 



iHQ 



N 



(32) 



As it is well known, C increases with the number of ex- 
cited channels Mi and the signal to noise ratio ~ 1/N 

M- 

In order to find the covariance matrix Q that maxi- 
mizes (I), we decompose the transmission matrix into a 
product of unitary matrices U and V and a M^ x Af 1 rect- 
angular diagonal matrix ^/T that contains the square- 
root of the transmission eigenvalues: 



t ^ uVfv. 



(33) 



Furthermore, we assume that the wave-guide is in the 
quasi- ID geometry L 3> VF, in which we can use the 
approximation 



dp{i) a p{{T,,})U„dT,A^i[U]d^i[V], 



(34) 



where d/i[t/] and d/x[y] are invariant measures of the 
unitary group. This approximation is known as the 
"isotropy assumption" since it implies that the flux in- 
cident in one scattering channel is, on av erag e, equally 
distributed among all outgoing channels [10[. On the 
other hand, since Q is positive definite, it can be writ- 
ten as Q = WDW\ where W is unitary and D diagonal 
and positive definite. Thus, the matrix product that ap- 
pears in Eq. jSOl) is tQP = uVt{VW)D{VW)1VtW. 
According to Eq. ^^, p{t) is unchanged under the trans- 
formation t — >■ tW, meaning that Q can be considered as 
a diagonal positive definite matrix without loss of gener- 
ality. We now adapt a proof given in Ref. 12| , initially 
proposed for Gaussian transmission matrices. First, we 
build the matrix 



Ml! 



(35) 



where {$«} is the ensemble of diagonal permutations ma- 
trices of size Ml X Ml. By construction, Q is equal to 
Tr(Q)/Mi/Mi. Then, we use the fact logjdet is concave 
on the set of positive definite matrices, so that the mu- 
tual information I{Q) considered as a function of Q is 
also concave. By virtue of Jensen's inequality, X{Q) ^ 



J2n^i^nQ^n)/Mil. Finally, using the isotropy approx- 
imation p4l) , we note that p{t ) is unchanged under the 
transformation i — >• i$„, so that (Z($„Q$„)) = (Z(Q)). 
Hence, {I{Q)) ^ {X{Q)) . We conclude that the matrix 
Q that maximizes {X[Q)) is proportional to the identity 
I Ml , which leads with (^5]) to 



Q = F'"7, 



Ml 



(36) 



Combining this result with Eq. (|32|) . we obtain Eq. (11) 
of the MT. 

In a waveguide with a geometry that does not satisfy 
L » W, the isotropy assumption ([M|) might still be valid 
in the bulk of the medium, if the wave propagation is far 
from the localized regime [lO|. However, the conversion 
of free propagation into diffusion at the boundaries of 
the medium makes that all channels do not have, on av- 
erage, the same transmission (see discussion above). It 
indicates that the distribution pit ) is not invariant under 
the permutation f — )• i <!>„, so that the proof of Eq. ([36)) is 
not entirely valid. That said, even in this situation, the 
covariance (p6| still leads to a good estimate of the typ- 
ical achievable capacity: it is the most reasonable choice 
when the statistical properties of t are not known at the 
receiver and, in fact, because of its simplicity and be- 
cause many space-time coding schemes conform to it, it 
is appealing even if some statistical information about t 
is available [2|. 
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